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X ■ 1 Introduction 

The operators of fractional differentiation have been introduced by Leib- 
nitz, Liouville, Riemann, Grunwal and Letnikov [6J. The fractional deriva- 
tives and integrals are used in the description of some models in mechanics, 
physics jGj, economics |1] and medicine [II]. The fractional variational cal- 
culus [1] is an important instrument in the analysis of such models. The 
Euler-Lagrange equations are non-autonomous fractional differential equa- 
tions in those models. 

In this paper we present the fractional jet fibre bundle of order A; on a dif- 
ferentiable manifold as being J"^(M, M) = Mx C'sc"'=(M), a E (0, 1), k G M*. 
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The fibre bundle 3°"^ is built in a similar way as the fibre bundle by R. 
Miron [S]. Among the geometrical structures defined on J°(R, M) we con- 
sider the dynamical fractional connection and the fractional Euler-Lagrange 
equations associated with a function defined on J"^(R, M). 

In section 2 we describe the fractional operators on M and some of their 
properties which are used in the paper. In section 3 we describe the frac- 
tional osculator bundle of order k. In section 4 the fractional jet fibre bundle 
J"(]R, M) is defined, the fractional dynamical connection is built and the 
fractional Euler-Lagrange equations are established using the notion of frac- 
tional extremal value and classical extremal value on J"'^(]R, M). In section 
5 we consider some examples and applications. 

2 Elements of fractional integration and dif- 
ferentiation on M 

Let / : [a, 6] ^ M be an integrable function and a G (0, 1). The left-sided 
(right-sided) fractional derivative of / is the function 



-Dtf)it) — r(l-a) dt la '^V-s')" ^ 

1 d fib)~f(s) . 



(1) 



y+^t J){T^) - r(i-a) dt it {s~tr 

where t G [a, h) and F is Euler's gamma function. 

Proposition 1. (see [6j) The operators ^Df and have the properties: 

1. If /i and /2 are defined on [a, 6] and exists, then 

-Dt{c,h + C2/2)(t) = Cl(-A"/l) W + C2(-A"/2) W- (2) 

2. If {an}n>Q is a real number sequence with lim = 1 then 

n— >oo 

lim(_D-/)(t) = (-AV)(t) = ^/(t)- (3) 

n^oo at 

3. a) If fit) = c,te [a, 6], c G M then 

(-A°/)W = o. (4) 

b) If f{t) = r',te [a, 6], 7 G M, then 
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c) If fit) = then 

= 1- (6) 
4. If fi and /2 are analytic functions on [a, b] then 

(-A"(/l/2))W=E| , I (-A"-Vl)Wpyi/2(t), (7) 



fc=0 



where (|r = |o|o...o|. 

5. It also holds true 

' h{t){^Dthmdt = - f h{t){+Dth){t)dt. (8) 

J a 

6. a) If / : [a, 6] ^ R admits fractional derivatives of order aa, a G N, then 

/(t+/i) = i?j(Mr„A")/w, (9) 

where Ea is the Mittag-Leffler function given by 

a=0 

b) If / : [a, 6] ^ M is analytic and G (a, 6) then the fractional McLaurin 
series is 



J-CtU 

fr^r(l + aa) 



i=0 
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The physical and geometrical interpretation of the fractional derivative 
on M is suggested by the interpretation of the Stieltjes integral, because 
the integral used in the definition of the fractional derivative is a Riemann- 
Stieltjes integral [TU] . 

By definition, the left-sided (right-sided) fractional derivative of /, of 
order a, m = [a] + 1, is the function 



D?m = ariL ^ds, G (-oo,t) 
*D?m = ^ i-irir^ds, G (t,oo). 



'121 



If supp/ C [a, b], then Dff = _Dff, *Dff = ^Df f . 
Let us consider the seminorms 



I tI — II * n'^T-ii 

I"''! ~ 11 '^IIl2(m) 1 



and the norms 

1/2 



1/2 



l-^lljgCM) ~ I ll^llL2( 



\x\ 


\J1{R)J 


\x\ 


lj«(R)) 



and J^i{W), J"^(]R) the closures of C(f^(M) with respect to the two norms 
from above, respectively. In [6J it is proved that the operators D° and 
satisfy the properties: 

Proposition 2. Let / C M and let J^i{I) and J^r{I) be the closures ofC^{I) 
with respect to the norms from above. For any x G J^^(/), < a < /3, the 
following relation holds: 

D^xit) = DtDf-xit). 
For any x G J^r{I), ^ < a < [3, it also holds 

*D^x{t) = *D^*Dfx{t). 
In the following we shall consider the fractional derivatives defined above. 



3 The fractional osculator bundle of order k 
on a different iable manifold 

Let a G (0, 1] be fixed and M a differentiable manifold of dimension n. 
Two curves p, a : I —>■ M., with p(0) = cr(0) = G M, G /, have a 
fractional contact a of order A; G N* in xq, if for any / G J^{U), xq E U, U a. 
chart on M, it holds 

Drifop)\t=o=Drifoa)\t=o (13) 

where a = 1, k. The relation f|T3l) is an equivalence relation. The equivalence 
class [p]°^ is called the fractional A;-osculator space of M in xq and it will be 
denoted by Osc'^^{M). If the curve p : I —>■ M is given by = t G /, 



class [p]°^ is called the fractional A;-osculator space of M in xq and it will be 

i = I, n, then, considering the formula ( fTTl) . the class [p]"^ , may be written 
as 



4.a 4.ak 

r(l + a) T[l + ak) 

where t G {—e,e). We shall use the notation 

x\Q)=x\ y^^'^-)= Drx\t)\t=o, (15) 

i (^i + aa) 
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for i = l,n and a = l,k. 

By definition, the fractional osculator bundle of order r is the fibre bundle 

{Osc''''{M),M) where Osc"^(M) = U^^^m Osc'^^{M) and 

Tvf : Osc'^''{M) ^ M is defined by 7r^''\[p]^^) = xq, (V)[p]^^*= G Osc°'=(M). 

For / G J^{U), the fractional derivative of order a, a E (0, 1), with 
respect to the variable x*, is defined by 

1 g rx^ /(xl,...,x'-l,s,a:'+l,...,x")-/(xl,...,a:'-l,a',x'+l,...,x") i ^ 
r(l-a) dx^ Ja^ (xi-s)" 

where are the coordinate functions on [/, z = 1, n, is the canonical 
base of the vector fields on U and Uab = {x E U, a* < < 6\ i = l,n} (Z U. 
Let [/, U' d M be two charts on M, [/ fl [/' 7^ and consider the change of 
variable 

= x*(x\...,x") (17) 

with det ^Ifj j 7^ 0. Let {c^a;*}^^— be the canonical base of 1-forms of V'^{U) 

and let us define the 1-forms (i(x*)" = a(x*)°~^(ix*, i = l,n. The exterior 
differential rf" : J^([/ n U') ^ V\U n U') is defined by 

(T = d{x^YD% = d{x^YD%. (18) 
Using (fT8l) and the property (^rjjif^j = 1, it follows that 

d{xr = T^77^,Du^rd{xr. m 

+ a) 

Using the notation 

Jli^, ^) = f(iT^^^-^'''^''' 

a 

d{x^Y = Jiix, x)d{xT- (21) 
From ([21]) it follows that 

^■(x, x)J^(x, x) = 5^ (22) 



from (]T9l) we get 



Consider x* = x*(t) and x*(t) = x*(x(t)), z = l,n, t G /. Applying the 
operator we get 

{D^x^m = D:,x\x){D^x^m = jp, x){D^x^m. (23) 
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Considering the notation from ffT^ we have 

y^H = Jj(x, (24) 

Also, from ([131) we deduce 

Dty'^-''^ = T^7^j^^,y'^'"'\ (25) 

1 [a[a — 1)) 



where i = l,n. Applying the operator in the relation fl24|) we find 

EM^yiiaa) ^ Y{1 + a) Jj(r^"~'^ x)?/^'(")+ 

riM ;j(y{"{-i)), y")2,.(2«) + ... + ;j(^«(-i), y"^)y.((^+iH+ (26) 

r(a(a-l)) ifaa) 

••• -T r(a) y ' 



where a = 1, A;. 

Proposition 3. (see |2], [S]) 

a) The coordinate transformation on Osc^"''\M), 

(x*, y*'^"'^^) — > (x*, ?/**^°\ 1/**^"^^) are given by the formulas ( fTTl) and 

(ESI). 



b) The operators D^i and the 1-forms ((ix*)", i = l,n, transform by the 
formulas 



a 



D^. = Jj(x, x)D:, ^^^^ 
dlx')" = Jj{x, x)d(xJ)". 

4 The fractional jet bundle of order k 

on a different iable manifold; geometrical 
objects 

By definition, the fc-order fractional jet bundle is the space J°''^(]R, M) = 
R X Osc''°'{M). A system of local coordinates on J"''(M, M) will be denoted 
by {t,x,y^°'\y^^"\ ...,y(^°)). Consider the projections : J"'=(R, M) ^ M 
defined by 

7r;^'=(t, X, y^"'^)) = X. (28) 
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Let U, U' C M be two charts on M with ?7n[/' ^ 0, (TT^y^U), K)"^(f/') C 
J"(M, M) the corresponding charts on J"(M, M) and, respectively, the cor- 
responding coordinates (x*), (x*) and (t, x\ y*^"^), {t, x\y^^"^). The transfor- 
mations of coordinates are given by 

/y»^ /y^ f rf^^ rft^^\ 

•AJ iX' I • • • • • ) 

a (29) 

Consider the functions (t)", (x*)", g J^((7r;^)~i(f/)), the 1-forms 

r(iT^«^(0", r(iT^«^(^0", r(iT^«^(?/*^"^)" ^ and the operators 



D", -D°i(„) on (ttq) i = l,n. The following relations hold: 

D'^.^^iw^M^^r) = TIT^d{t-)iD?) = 1, (30) 
On J"(R, M) we may define the canonical structures 



s = e^ 

a 

V- — D° 



(31) 



Using ( !29l) it is easy to show that the structures ( |3T1) have geometrical charac- 
ter. The space of the operators generated by the operators {Df, D'^i, D'^i^^-^}, 

i = 1, n, will be denoted by x"(('^o )~^(^))- ^'^^ a — > 1 the space of these 
operators represents the space of the vector fields on 71^^(17). 

a 

A vector field T G x"((7ro)~ (f^)) is called FODE (fractional ordinary 
differential equation) iff 

d{tY{T) = 1 

" ; (32) 

e\T) = 0, 



for i = l,n. In local coordinates FODE is given by 



T = D^ + + , (33) 
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where F' G C°°((7r(^)~^(t/)), i = l,n. The integral curves of the field FODE 
are the solutions of the fractional differential equation (EDF) 



Df"x'{t) = F\t,x{t),D^x{t)), i = l,n. (34) 

The fractional dynamical connection on J"(M, M) is defined by the frac- 
tional tensor fields H of type (1,1) which satisfy the conditions 

a a 

eioH = o 

e2oH = e2 (35) 

a 

H 



c = —id 

V 



V 



where V is formed by operators generated by {-D^,{c)}j=T7l- the chart 
(ttq )~ (t/) the fractional tensor field H has the expression 

1 2 . 3 ^ 

H = {H d{tY + H] d{x'Y + Hi d(y*("))°) ® 

4 5 6 

{Hi{dtY + H] d^x^Y + H] d(y^("))") ® D° + (36) 

7 8 9 

{w d{tY + H] d{x^Y + m c/(2/^("))") ® ) . 

J J y 

a 

The tensor field H has a geometrical character, fact which results by using 
the relations (12^ . and is called a (i°-tensor field. Using the relations 
and (pTl) we get 



Proposition 4. a) The fractional dynamical connection H, in the chart 
{'Kq)^^{U), is given by 

H = f(iT^[(-2/^^")^." + HW-,^,) ® d{tr+ ^^^^ 
{D^. + (.)) ® d{xY - ® d(|/^("))-]. 

a 

h) The fractional dynamical connection H defines a /(3, — 1) fractional struc- 
ture on J"(M, M), I.e., [h^ = H- 

a Q, 

c) The fractional tensor fields I and m which are defined by 

a a a 

l=HoH 

a a a (38) 

m = -HoH+L 
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where I is the identity map, satisfy the relations 

a a a ce a a c c a 

lol^l, mom^mol, i+m^I 

= /(^,".(.)) = ^;.(.) ^ (39) 

a 

c?j The fractional vector field V £ x"(J"(]R, M)) given by 

f = m( A") = A" + y'^"^^:^ + (2/^^"^ + ^'■)^,° (.) (40) 

defines a field FODE associated to the fractional dynamical connection. The 
integral curves are the solutions of the EDF 

Dl^x^t) = D'^x'it) Hi +r(l + a) H^ (41) 

where Hi and W are functions of {t., x{t), y^"'\t)). 

Let L e C°°( J"(R, M)) be a fractional Lagrange function. By definition, 

a 

the Cartan fractional 1-form is the 1-form 9l given by 

eL = Ld{t)'' + S{L). (42) 

We call the Cartan fractional 2-form, the 2- form ujl given by 

(^L = c?" Ol (43) 
where (i" is the fractional exterior differential: 

= (i(0" A" + ^^(^T^:* + d{y'^''YD^i(c) ■ (44) 
In the chart (tTq)" (t/), 6*^;, and uj^ are given by 

COL = Ac^(t)" A d(xO° + Bidif") A (i(i/*(°))"+ (45) 
Ajjd(a;*)" A (/(a;^)" + Bijd{xY A d(y^("))", 
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where 

Proposition 5. If L is regular (i.e., det (^ g^i(ffj^j(c.) j 7^ ^^en i/iere exzsfo 

a fractional field FODE Tl such that i a ujl = 0. In the chart {tcq)^ (U) we 
have 

Tl = Dt + + , (47) 

where 

df = Df + (48) 
(^7^'=) = (D-(„)DV,(L))-^ 

An important structure on J"(M, M) is described by tlie fractional Euler- 
Lagrange equations. Let c : t E [0,1]^ (^*(^)) € Af be a parameterized 
curve, such that Imc G U G M. The extension of the curve c to J"(M, M) 
is the curve : t e [0,1] (t, x*(t), y^(")(t)) G J"(M, M). Consider 
L G C°°(J"(M, M)). The action of L along the curve c'^ is defined by 

A{c") = f\{t,x{t),y"{t))dt. (49) 
Jo 

Let Cs : t E [0, 1] {x^{t, e)) G Af be a family of curves, where e is sufficiently 
small so that ImCe C U, Co{t) = c{t), DfCe{0) = DfCe{l) = 0. The action of 
L along the curves q is 

A{c^)= f L{t,x{t,e),y''{t,e))dt, (50) 

^0 

where y^^^^t, e) = p^^:^Df x*(t, e). The action (J50l) has a fractional extremal 
value if 

D^Aic^) |,=o = 0. (51) 



The action (l50l) has an extremal value if 

DlAic'^) |,=o = 0. (52) 
Using the properties of the fractional derivative we obtain 
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Proposition 6. a) A necessary condition for the action flHIj) to reach a 
fractional extremal value is that c{t) satisfies the fractional Euler-Lagrange 
equations 



(53) 



yi.{a) -I 



where i = l,n. 

b) A necessary condition for the action f lHIl) to reach an extremal value is 
that c{t) satisfies the Euler-Lagrange equations 



(54) 



,,i(a) 1 



where i = l,n. 

The equations fl53p may be written in the form 

D^.L - dtiD^.,^,L) - 2/^(2°) D;,,., (Z^;.(.)L) = 0, (55) 



for i = l,n. The equations (ImII may be written as 



where z = 1, n. Let us denote by 

9g = D°,.,(DJ,.,L), (57) 

and by = (,g)-', if det(,g) ^ 0. Fro.n 63 and from ProposMon 5, 

a 

we get the fractional field FODE associated to L. 

Let c : t e [0, 1] — (x*(t)) C ?7 be a parameterized curve. The extension of 
c to J^(M, M) is the curve c"'= : t G [0, 1] ^ (t, x*(t), ?/""(t)) G J"*=(M, M), 
a = 1, /c. Let L : J°'^(M, M) — M be a Lagrange function. The action of L 
along the curve c°'^ is 

A{c"'')= [ L{t,x{t),y''\t))dt. (58) 

Let Ce : t E [0,1] — > (a;*(t,£)) G M be a family of curves, where the abso- 
lute value of e is sufficiently small so that ImCe C U C M, co(t) = c(t), 
D"c{e) |e=o = -D" c(£) |e=i = 0. The action of L on the curve q is given by 



A{cf)= l\{t,x{t,e),y"%t,e))dt (59) 
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where y^^""^\t, e) = pp-|^Z^f"x*(t, e), a = l,k. The action fl3U|) has a frac- 
tional extremal value if 

D^iAicf)) |,=o = 0. (60) 

The action fl59l) has an extremal value if 

Dl{A{cf)) |,=o = 0. (61) 

Proposition 7. a) A necessary condition for the action fl3H|) to reach a 
fractional extremal value is that c{t) satisfies the fractional Euler-Lagrange 
equations 

k 

d:.l + Y,i-^rdr{D^,^..^L) = o, (62) 

a=l 

where 

dr = + + 2/^(2") + ... + y^^^'^^D" (63) 



and i = l,n. 

b) A necessary condition that the action ( l58l) reaches an extremal value is 
that c{t) satisfies the Euler-Lagrange equations 

— + 5^(-l)X(^,V)^) = 0, (64) 

a=l 

where 

d^ = D] + y^^^^Dl. + ... + (65) 
Example. Consider the fractional differential equation 



a2r(l + 3a)y(=^") = 0. 



(66) 



The equation ( 166|) is the fractional Euler-Lagrange equation ( l62l) for the 
function 

L = j^x'^-a,T{l + 2a){y^r+ 
a2r(l + 3a)(2/2")-. 

The equation (lUUj) is the fractional Euler-Lagrange equation for the 
function 

L = + / - -r(i + 2«)(y")^ + ^r(i + 3a)(,-)^ 
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5 Examples and applications 

1. The nonhomogeneous Bagley-Torvik equation 

The dynamics of a flat rigid body embedded in a Newton fluid is described 
by the equation 

aD^x{t) + hDl'^x{t) + cx{t) - f{t) = 0, (67) 

where a, 6, c G M and the initial conditions are x{0) = 0, Dlx{0) = 0. The 
equation flU7|) is a fractional differential equation on the bundle J°(R, M) for 
ex = \. Indeed, let's consider the fractional differential equation 

aDl'^xit) + bDl'^xit) + exit) - f{t) = 0, (68) 

with a > 0. For a = ^ the equation ( l68l) reduces to (1671) . With the notations 
(fT5|) . the equation (!68l) becomes 



ar(l + 8a)?/(^")(t) + 6r(l + 6a)2/(^")(t) + cx(t) - /(t) = 0. (69) 
On the bundle J'^°'(]R, M) let us consider the Lagrange function 



L(t,x,y(3")^^(4«)) = icx2 -/x- |r(l + 6a)(2/(3"))2+ 
fr(l + 8a)(2/(4' 



q)\2 



(70) 



Using the relation (1651) . the Euler-Lagrange equation for ( ITOll is 

cx- f + bT{l + 6a)D3"?/(3") + ar(l + 8a)D^'^y^^''^ = (71) 
cx - / + br(l + 6a)?/(6°) + aT{l + 8a)y^^"^ = 0. 



Proposition 8. The equation ( 1671) represents the Euler-Lagrange equation 
on the bundle J^"(]R, R) for a=\, with the Lagrange function given by 

L{t,x,y(y'\y(')) = \cx^ - fx-\T{^/2){y^^l^)f+ ^^^^ 
fr(3)(y(2))2. 

2. Differential equations of order one, two and three which admit 
fractional Lagrangians 

The following differential equations don't have classical Lagrangians such 
that the Euler-Lagrange equation represents the given equation: 

x(t) + Vi(t,x) = 0, Vr{t,x)= \ (73) 
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x{t) + a,xit) + V2it,x) = 0, V2it,x) = ^^f'''\ (74) 

ox 

x{t) + a2x{t) + aix{t) + V3{t,x) = 0, V3{t,x) = (75) 



Let us associate the fractional equations from below to the equations ( !73l) . 
( |7il) and ( 1751) . respectively: 

D2"a;(f) + y^(f,a;) = 0, (76) 

Df'^x{t) + aiDl^'xit) + V^2(t, x) = 0, (77) 
A^°x(t) + a2Df''x{t) + aiDl^xit) + V^it, x) = 0. (78) 

Proposition 9. a) Let J"(M, M) ^ M &e the fractional bundle and consider 
L : J°(M, M) ^ M ^zven by 

L(t,x, ?/(")) = t/i(t,x) - ^r(l + 2«)(2/")^ (79) 
r/ie Euler- Lagrange equation of ( I79l) zs 

S-A"(^) = ^ + r(i + 2a),(-) = 



V^i(t,x) + D2"2;(t) =0. 



^0) 



6^ Le^ J^"(M, M) — > M &e the fractional bundle and the Lagrangian L : 
J2"(M, M) M gzven by 



L{t, X, ?/(2a)) _ ^) _ + 2a)(|/")2 + 

ir(l + 4a)(y(2"))2. 
r/ie Euler- Lagrange equation of ( IHTi) 



^1) 



^2(^,0;) +air(l + 2a)?/(2")+ 

a2r(l + 4a)2/(^") = 

V2{t,x) + aiD2°x(t) + A^°x(t) = 0. 



^2) 



Let J3°(]R, M) ^ M 6e i/ie fractional bundle and L : J3"(M, M) ^ M given 

by 

L(t,x, ?/(2")^^(3a)) = V3(t,x) - f r(l + 2a)(?/(°))2+ 
f r(l + 4a)(y(2"))2 _ ir(i + 6a)(y(3"))2. 



83) 
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The Euler- Lagrange equation of flH^ is 

f - D? (^) + (^) - Dt (^) = V,it, .) + 

aiT{l + 2a)y(2") + ^21(1 + 4a)?/(4") + r(l + 6a)y^^'^'> = (84) 

V3{t,x) + aiL>2a^(t) + a2£'f"x(t) + = 0. 

d) For a = ^ we obtain the fractional Lagrangians that describe the equations 
(I7H|) . dH]), (USD, respectively 

L{t,x,yW^^) = U^{t,x) - '^T{2){yW^)Y 

L(t,x,y(V2),^(i)) = - iair(2)(y(V2))2 + lr(3)(y«)2 

L(t,x,y(V2), ^(1)^^(3/2)) ^ _ f r(2)(2/(V2))2+ 

fr(2)(2/«)2-ir(4)(y(3/2))2. 

In the category of the equations (17^ and (17^ there are: 

a) the nonhomogeneous classical friction equation 

mx(t) + jx(t) - = 0, (86) 

ox 

b) the nonhomogeneous model of Phillips [8] 

x(t) + aii;(t) + 6ix(t) + /(t) = 0, (87) 

c) the nonhomogeneous business cycle with innovation [8] 

y{t) + a2y{t) + aiy{t) + hx{t) + /(t) = 0. (88) 

Conclusions 

The paper presents the main different iable structures on J"(M, M), in order 
to describe fractional differential equations and ordinary differential equa- 
tions, using Lagrange functions defined on J"(M, M). 

With the help of the methods shown, there may be analyzed other models, 
such as those found in [1] and [TT] . 
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